Immune-mediated competition in rodent malaria is most likely caused by induced changes in innate immune clearance of merozoites

Malarial infections are often genetically diverse, leading to competitive interactions between parasites. A quantitative understanding of the competition between strains is essential to understand a wide range of issues, including the evolution of virulence and drug resistance. In this study, we use dynamical-model based Bayesian inference to investigate the cause of competitive suppression of an avirulent clone of Plasmodium chabaudi (AS) by a virulent clone (AJ) in immuno-deficient and competent mice. We test whether competitive suppression is caused by clone-specific differences in one or more of the following processes: adaptive immune clearance of merozoites and parasitised red blood cells (RBCs), background loss of merozoites and parasitised RBCs, RBC age preference, RBC infection rate, burst size, and within-RBC interference. These processes were parameterised in dynamical mathematical models and fitted to experimental data. We found that just one parameter μ, the ratio of background loss rate of merozoites to invasion rate of mature RBCs, needed to be clone-specific to predict the data. Interestingly, μ was found to be the same for both clones in single-clone infections, but different between the clones in mixed infections. The size of this difference was largest in immuno-competent mice and smallest in immuno-deficient mice. This explains why competitive suppression was alleviated in immuno-deficient mice. We found that competitive suppression acts early in infection, even before the day of peak parasitaemia. These results lead us to argue that the innate immune response clearing merozoites is the most likely, but not necessarily the only, mediator of competitive interactions between virulent and avirulent clones. Moreover, in mixed infections we predict there to be an interaction between the clones and the innate immune response which induces changes in the strength of its clearance of merozoites. What this interaction is unknown, but future refinement of the model, challenged with other datasets, may lead to its discovery

Ecological divergence and hybridization of neotropical leishmania parasites

The tropical Andes are an important natural laboratory to understand speciation in many taxa. Here we examined the evolutionary history of parasites of the Leishmania braziliensis species complex based on whole-genome sequencing of 67 isolates from 47 localities in Peru. We first show the origin of Andean Leishmania as a clade of near-clonal lineages that diverged from admixed Amazonian ancestors, accompanied by a significant reduction in genome diversity and large structural variations implicated in host-parasite interactions. Within the Andean species, patterns of population structure were strongly associated with biogeographical origin. Molecular clock and ecological niche modeling suggested that the history of diversification of the Andean lineages is limited to the Late Pleistocene and intimately associated with habitat contractions driven by climate change. These results suggest that changes in forestation over the past 150,000 y have influenced speciation and diversity of these Neotropical parasites. Second, genome-scale analyses provided evidence of meiotic-like recombination between Andean and Amazonian Leishmania species, resulting in full-genome hybrids. The mitochondrial genome of these hybrids consisted of homogeneous uniparental maxicircles, but minicircles originated from both parental species. We further show that mitochondrial minicircles-but not maxicircles-show a similar evolutionary pattern to the nuclear genome, suggesting that compatibility between nuclear-encoded mitochondrial genes and minicircle-encoded guide RNA genes is essential to maintain efficient respiration. By comparing full nuclear and mitochondrial genome ancestries, our data expand our appreciation on the genetic consequences of diversification and hybridization in parasitic protozoa.status: publishe

A theoretical study of random segregation of minicircles in trypanosomatids

The kinetoplast (k) DNA network of trypanosomatids is made up of approximately 50 maxicircles and the order of 10(4) minicircles. It has been proposed, based on various observations and experiments, that the minicircles are randomly segregated between daughter cells when the parent cell divides. In this paper, this random segregation hypothesis is theoretically tested in a population dynamics model to see if it can account for the observed phenomena. The hypothesis is shown to successfully explain, in Leishmania tarentolae, the observation that there are a few major and many minor minicircle classes, the fluctuations of minicircle class copy numbers over time, the loss of non-essential minicircle classes, the long survival times of a few of these classes and that these classes are likely to be the major classes within the population. Implications of the model are examined for trypanosomatids in general, leading to several predictions. The model predicts variation in network size within a population, variation in the average network size and large-scale changes in class copy number over long time-scales, an evolutionary pressure towards larger network sizes, the selective advantage of non-random over random segregation, very strong selection for the amplified class in Crithidia fasciculata if its minicircles undergo random segregation and that Trypanosoma brucei may use sexual reproduction to maintain its viability

Organization of minicircle cassettes and guide RNA genes of Trypanosoma brucei

Mitochondrial DNA of protists of order Kinetoplastida comprises thousands of interlinked circular molecules arranged in a network. There are two types of molecules called minicircles and maxicircles. Minicircles encode guide RNA (gRNA) genes whose transcripts mediate post-transcriptional editing of maxicircle encoded genes. Minicircles are diverse. The human sleeping sickness parasite Trypanosoma brucei has one of the most diverse sets of minicircle classes of all studied trypanosomatids with hundreds of different classes, each encoding one to four genes mainly within cassettes framed by 18 bp inverted repeats. A third of cassettes have no identifiable gRNA genes even though their sequence structures are similar to cassettes with identifiable genes. Only recently have almost all minicircle classes for some subspecies and isolates of T. brucei been sequenced and annotated with corresponding verification of gRNA expression by small-RNA transcriptome data. These data sets provide a rich resource for understanding the structure of minicircle classes, cassettes and gRNA genes and their transcription. Here, we provide a statistical description of the functionality, expression status, structure and sequence of gRNA genes in a differentiation-competent, laboratory-adapted strain of T. brucei. We obtain a clearer definition of what is a gRNA gene. Our analysis supports the idea that many, if not all, cassettes without an identifiable gRNA gene contain decaying remnants of once functional gRNA genes. Finally, we report several new, unexplained discoveries such as the association between cassette position on the minicircle and gene expression and functionality, and the association between gene initiation sequence and anchor position

Estimates of , in reconstituted (6 mice), nude (7 mice) and wildtype (6 mice) in single and mixed infections. Bars represent the means of the means of the marginal posteriors.

<p>Error bars represent 1 standard error.</p

Model fits to parasite densities.

<p>Fits of the single-cause model (H₅) with clone-specific to AS (green) and AJ (blue) parasite densities in reconstituted (top panels), nude (middle panels) and wildtype (bottom panels) mice during mixed infections. Crosses are data. The solid lines give the median fits. Grey regions correspond to the 95% posterior intervals of model uncertainty. These plots show that the model fits the data quite well for each individual in all treatment groups.</p

Assessment of functional forms of the adaptive immune responses.

1<p>Mean2 standard errors of five independent fits.</p